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Abstract 

How the vortex lattice orders at long range in a layered superconductor with weak point pinning 
centers is studied through a duality analysis of the corresponding frustrated XY model. Vortex- 
glass order emerges out of the vortex liquid across a macroscopic number of weakly coupled layers 
in perpendicular magnetic field as the system cools down. Further, the naive magnetic-field scale 
determined by the Josephson coupling between adjacent layers is found to serve as an upperbound 
for the stability of any possible conventional vortex lattice phase at low temperature in the extreme 
type-II limit. 



INTRODUCTION 



It is well known that an external magnetic field can penetrate a type-II superconductors in 
the form of lines of flux quanta 1]. The repulsive forces that such flux lines experience favor 
the creation of a triangular vortex lattice, while the quenched point disorder present to some 
degree in all superconductors frustrates that tendency. Three thermodynamic groundstates 
are then likely. Either the triangular vortex lattice is robust to weakpoint pinning and 
assumes a Bragg glass state with no lines of dislocations that thread it|2j, or it will transit 
into a defective state with quenched-in lines of dislocations that thread it. The latter, in turn, 
has two possible outcomes: a vortex glass state that retains macroscopic phase coherence of 
the superconducting order parameter^] , or a pinned liquid state that does notQ]. 

High-temperature superconductors, in particular, are extremely type-II and layered^]. 
Below, we shall study how a vortex lattice pinned by material point defects orders at long- 
range in such materials. The vortex lattice in layered superconductors with weak random 
point pins shall be described theoretically in terms of the phase of the superconducting order 
parameter via the corresponding frustrated XY model 4]. This model notably neglects the 
effects of magnetic coupling between layers, while it treats the Josephson coupling between 
them exactly. The growth of long-range order across layers is then computed from the XY 
model through a duality analysis[5], where the ratio of the energy of the Josephson coupling 
between adjacent layers to the temperature emerges as a small parameter. We find first 
that the correlation length for vortex-glass order [^J across weakly coupled layers diverges as 
temperature cools down from the vortex liquid towards the two-dimensional (2D) ordering 
transition. The divergence signals a transition to a vortex glass phase jf]]!?]]!^^. Second, 
we find no evidence for the divergence of conventional superconducting phase correlations 
across layers from inside the latter vortex glass to lowest order in the inter-layer Josephson 
coupling. This indicates ultimately that the naive decoupling field for the pristine vortex 
lattice serves as an upper bound for a stable Bragg glass phase |2j in the extreme type-II 
limit. Comparisons with previous numerical [J], theoretical U ], and experiment alQ] deter- 
minations of the stability line for the Bragg glass in layered superconductors are made at 
the end of the paper. 
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TWO DIMENSIONS 



The XY model with uniform frustration is the minimum theoretical description of vortex 
matter in extremely type-II superconductors. Both fluctuations of the magnetic induction 
and of the magnitude of the superconducting order parameter are neglected within this 
approximation. The model hence is valid deep inside the interior of the mixed phase. The 
thermodynamics of an isolated layer with uniform frustration is determined by its superfluid 
kinetic energy 

e xy = - E E -W a m* - M* C 1 ) 

which is a functional of the phase of the superconducting order parameter, e^, over the 
square lattice, r. Here, J x and J y are the local phase rigidities that are equal and constant, 
except over links in the vicinity of a pinning center. The vector potential A = (0,2irfx/a) 
represents the magnetic induction oriented perpendicular to the layers, B± = Q f/a 2 . Here 
a denotes the square lattice constant, which is of order the coherence length of the Cooper 
pairs, $o denotes the flux quantum, and / denotes the concentration of vortices per site. 

Analytical and numerical work indicates that the 2D vortex lattice is invaded by 
quenched-in dislocations in the presence of any degree of random point pinning jl^] . The 
author has argued [l^ that the dislocations quenched into each 2D vortex lattice described 
by the frustrated XY model (PJ notably do not line up to form low-angle grain bound- 
aries, however (cf. ref. |l5(]). That argument is based on the incompressible nature of 2D 
vortex matter in the extreme typ e-II limit. The absence of grain boundaries is consistent 
with Monte Carlo simulations |l6( of the equivalent 2D Coulomb gas ensemble with random 
point pins[8j], as well as with Monte Carlo simulations of the frustrated XY model in three 
dimensions with randomly located columnar pins 17]. Secondly, a net superfluid density is 
predicted at zero temperature for perpendicular magnetic fields above the collective-pinning 
threshold, B^ 2D \ in which case the number of pinned vortices is greater than the num- 
ber of isolated dislocations quenched into the 2D vortex lattice 18|. Here, the scale of the 
Larkin domains is set by the separation between neighboring dislocations quenched into 
the vortex lattice. A variational calculation by Mullock and Evetts yields the estimate 
B^' ~ (4f p /Eod) 2 &o for the threshold field where f p denotes the maximum pinning 
force, where 6q = (^o/AtcXl) 2 is the maximum tension of a fluxline in the superconductor, 
and where d denotes the separation between adjacent layers. Here Al represents the London 
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penetration depth. The pinning of the vortex lattice in isolated layers shall be assumed to 
be collective henceforth: B± > Bj? D \ 

The previous indicates that a hexatic vortex glass characterized by a homogeneous dis- 
tribution of quenched-in dislocations and by a net superfluid density exists in isolated layers 
of the frustrated XY model with weak random point pins at zero temperature |l8j| . 
The transition temperature T^ 2D ^ that separates the low-temperature hexatic vortex glass 
from the high-temperature vortex liquid must therefore be equal to zero or greater. Recent 
current-voltage measurements of 2D arrays of Josephson junctions in weak external mag- 
netic field indicate that the 2D superconducting/normal transition at T = T^ 2D ^ is second 
order [2^, with Tj 2D ) much larger than the 2D melting temperature of the pristine vortex 
lattice, T^ D ) = J/20. Since the previous is a faithful realization of the frustrated XY model 
in 2D with random point pinning centers, we shall assume henceforth that the hexatic 
vortex glass melts into a vortex liquid at temperature Tj 2D ) > via a second-order phase 
transition. 



THREE DIMENSIONS 



We shall now demonstrate how long-range vortex-glass order emerges across layers from 
the vortex liquid phase of layered superconductors with weak random point pins. Let us first 
couple the layers through the Josephson effect by adding a term — J z cos(A 2 — A z ) to the 
internal energy of the frustrated XY model for each nearest-neighbor link across adjacent 
layers. The component of the magnetic induction parallel to the layers is taken to be null 
throughout. At weak coupling, J z <C IcbT, phase correlations across N layers can then be 
determined from the quotient 

(exp[i5>(r)0(r)]) = Z CG [p] /Z CG [0] (2) 

r 

of partition functions for a layered Coulomb gas (CG) ensemble fl: 

Zcg[ P ]= E yo^-^iQlq^e-^r^. (3) 

Above, n z (f,l) is a dual charge/integer field that lives on links between adjacent layers I 
and I + 1, located at 2D points r, and p(r) = <5 r -? • (<^,i — S^n) is the external integer probe 
field. The ensemble is weighted by a product of phase auto-correlation functions for isolated 



layers I, 

Q[q] = (exp[i£g(fW, Q]>j,=o, (4) 
probed at the dual charge that accumulates onto that layer: 

Qi(r) = p(r, I) + n z (r, I - I) - n z (r, I). (5) 

It is also weighted by a bare fugacity yo that is raised to the power iV[n 2 ] equal to the total 
number of dual charges, n z = ±1. The fugacity of the dual CG ensemble is given by 
Ho = Jz/SksT in the selective high-temperature regime, J z <C ksT, reached at large model 
anisotropy. It is small compared to unity in such case, which implies a dilute concentration 
of dual n z charges [5]. The dual CG ensemble (j3J) is valid in that regime. 

The above duality analysis is particularly natural and effective in the vortex-liquid phase, 
where autocorrelations of the superconducting order parameter in isolated layers (J1J) are short 
range. They shall be assumed to take to the form that is characteristic of a hexatic vortex 
liquid between points f\ and f 2 in an isolated layer /[3]j2l|: 

Ci(l, 2) = g oe - r ^/^e- iMl) e iM2) . (6) 

Here e l< ^° is the superconducting order parameter of layer I in isolation at zero temperature, 
^2D denotes the phase correlation length of the 2D hexatic vortex liquid, and go is a prefactor 
of order unity. Also, r 12 = r i — r : 2 is the displacement between the probes within layer I. 
To lowest order in the (dual) fugacity, yo, Eqs. (J2J) and © then yield the expression 



(e#M + n) = y^-J2 Ci(0, 1) • C m (l, 2) • ... • C l+n (n, 0) (7) 

1 n 

for the bulk average (overbar) of the gauge-invariant auto-correlation function of the con- 
ventional superconducting order parameter e l< ^ across n layers, at zero parallel field. Above 
and hereafter, we take the gauge A z = 0. The uncorrelated nature of point pinning centers 
across layers implies the form 



]J^_ e~ i ^o(r m ,l+rn) e i<t> (r m+1 ,l+m) _ n^ l _ e _rm,m+1 ' * (8) 

for the bulk average of the relevant product of zero-temperature order parameters, with 
matching endpoints r = r n+1 . Here, 1$ is a quenched disorder scale that is set by the 
density of lines of dislocations quenched into the vortex lattice at J z = that begin or end 
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at a given layer. We remind the reader that 1$ is believed to be finite (in the absence of 



inter-layer coupling) for any non-zero strength of quenched point disorder [13|. Substitution 
of (JHJ) into expression (J7J) then yields the principal dependence|22] 



(e^M + n) oc [g (J/k B T)((l^ + tfr^Ao) 2 )" (9) 

for the correlation of the conventional superconducting order parameter across n layers. 
Here, J is the macroscopic phase rigidity of an isolated layers at zero temperature, Ao = 
(J/Jz) 1 ^ a is the Josephson penetration depth, and ^ = ^2d/2. Notice that the existence 
of the disorder scale 1$ implies that the perturbative result Q above does not diverge with 
the 2D phase correlation length in the vicinity of the 2D ordered phase. We conclude 
that conventional superconducting phase coherence across many layers (n — > oo) does not 
emerge out of the vortex liquid at weak Josephson coupling between adjacent layers. 

The growth of macroscopic vortex-glass order across layers Q from inside of the vortex 
liquid is still possible, however. We shall test for it by computing the corresponding auto- 
correlation function |3|, which is given by 



|( e #M + «)|2 - y g» J2 ... £ Ci(0, 1)Q(0, 1) • a+i(l, 2)Q +1 (1, 2) ■ ... ■ a+n(n, 0)Q +n (n, 0) 

1,1 n,n 

(10) 

to lowest order in the (dual) fugacity, Uq. It is natural to look for vortex-glass order to 
emerge from within the 2D critical regime: £ 2D > 2Z at T > T^ 2D \ where denotes 
the transition temperature of the 2D hexatic vortex glass. The bulk average of the product 
of zero-temperature order parameters that appears in the integrand above can then be 
approximated by the corresponding product of the bulk averages limited to adjacent layers, 
V — I + m — 1 and V + 1, only: 



exp[#{,°J, +1 (m)] • expH4°], +1 (m)] = e~ r ^^. (11) 

Here = 4>o(r, I' + 1) — 0o(^*? — A z (r) is the quenched inter-layer phase difference. 

Converting to center-of-mass variables among the inter-layer coordinates, f m and f m , then 
yields the principal dependence!^ 



|<e^+">| 2 ex \go(J/k B T){l&/l4)r ( 12 ) 

for the vortex-glass correlations across layers in the 2D critical regime, £20 S> 2/^, at zero 
parallel field. The corresponding correlation length is equal to the layer spacing (d) 



when the argument in brackets above is set to 1/e. This occurs at a cross-over field 

B x ~ goiJ/ksm^/aD^o/Al) (13) 

that separates two-dimensional from three-dimensional (3D) vortex-liquid behavior (see Ta- 
ble H}. Above, a vx denotes the square root of the area per vortex inside of a given layer. Also, 
the argument between brackets on the right-hand side of Eq. (|12p notably diverges with 
in the vicinity of the 2D ordering transition. This indicates that a transition to a vortex 
glass that orders across a macroscopic number of layers [(|, £j_ — > oo, occurs at a critical 
temperature T g that lies inside of the window [T^ 2D \T X ]. Indeed, setting the argument of 
the exponent on the right-hand side of Eq. (|T2)l to unity yields a critical field B g = B x /e, 
below which a vortex glass exists (see Table [I]). 

Last, recall that the superfluid density across layers, = —Af~ 1 k B Td 2 \nZ C G/dA 2 .\ , is 
given by the expression 

^•^([E^OD^T, (14) 

r,l 

where M counts the number of nearest-neighbor links between layers, and where periodic 
boundary conditions are assumed across layers. Study of Eqs. ©-© yields that the tension 
for a line across layers of dual n z quanta is equal to £j x , where £j_ denotes the correlation 
length for vortex-glass order across layers. The corresponding superfluid density (pPfj) is then 
null in the limit of a macroscopic number of layers inside of the vortex liquid, where £j_ < oo 
(see Table IU). 

The previous result (J12|) clearly demonstrates that a selective high-temperature expansion 
in powers of the fugacity yo necessarily breaks down in the 2D ordered phase, T < T^ 2D \ 
where £20 is infinite. A direct analysis of the frustrated XY model for an isolated layer 
finds, in particular, that long-range correlations of the superconducting order parameter (0J) 



[8] 21]: 



decay algebraicly instead at such low temperatures 

Q[q] = g'o + ■ exp r) 2D J2 q(ri)\n(r l!2 /r )q(r 2 ) ■ exp[i J2 <?(ji)0o(n, I) . (15) 

(1.2) 1 

The exponent r] 2 D that characterizes the algebraic decay of 2D phase coherence is related 
to the 2D superfluid density by pf D ^ = k B T/2irri 2 D- Above, g = pf^/J is the ratio of the 
2D phase stiffness with its value at zero temperature, J, while n + counts half the number of 
probes in q(r). Also, r denotes the natural ultraviolet scale. It is important to observe at this 



stage that the loop excitations in the (completely) dual representation of the 3D XY model[5( 
lose their integrity in the ordered phase. This translates into the absence of charge conser- 
vation in the (partially) dual CG ensemble (jSJ). In other words, the dual n z charges form a 
plasma in the ordered phase. A Hubbard- Stratonovich transformation of the CG partition 
function (j3J) followed by the unrestricted summation of configurations of charges with values 
n z = 0, ±1 then yields the equivalent partition function 2^] Zld[p] = J V8 e^ EhD ^ kBT+l ^ p ' e 
for a renormalized Lawrence- Doniach (LD) model that shows no explicit dependence on the 
perpendicular magnetic field. Its energy functional is specifically given byp g| 



^(W,) 2 -Ao 2 cosfl M+1 



(16) 



where 6ij+i = 4>u+\ + ~ The above continuum description is understood to have an 
ultraviolet cut off r of order the inter- vortex spacing a vx . 

We can now determine the growth of correlations across layers of the conventional su- 
perconducting order parameter deep inside of the vortex glass phase, T < T^ 2D \ at weak 
Josephson coupling between layers, A —>■ oo. The physics described by the original layered 
XY model coincides directly with that of the renormalized LD model described above at 
large scales in distance compared to the ultraviolet cutoff, tq. Asymptotic correlations of 
the conventional superconducting order parameter across layers, for example, are identical 



to those of the LD model: lim rwoo (e l<? ^ ! + n ) = (e ! i+n ). The configuration that optimizes 
.Eld must be determined first in order to compute the later near zero temperature. The LD 
energy functional (fTtjj) implies that it satisfies the field equation 

- Vyfl - 6?] = A - 2 sm^ 1/+2 - 2Ar 2 sin^°i, +1 + A 2 sin ft,,, (17) 

where Ai = A (cf. refs. j^J and 2^)- The phase angles 9^ are then constant inside of a 
given layer /' in the weak coupling limitj^, A , A x — > oo. Next, if 58^ denotes the fluctuation 
in the phase angles, the auto-correlation function for conventional superconducting order 
across many layers is then approximated by the expression 



n(0) 



e ie u+n c* njtr e ■ w+i - m>] (is) 



near zero temperature, to lowest order in the fluctuation 26] . After inverting the field equa- 
tion (JT7j) for the fluctuation of the phase difference between adjacent layers, substitution 



S 



into the expression above yields the result 



a n A 1 



2n 



n^ =1 / d 2 r m G^(0,m) 



Yl^_ ie M+m-l,l+m(°) e r i ^i+m-l,l+m( m ') (19) 



for the autocorrelation of the superconducting order parameter across layers. The prefactor 
on the right-hand side satisfies the recursion relation a n+ i = a n + {K\/2K^) 2 a n _i, with a® = 1 
and a_i = 0. Also, 

G m = [_ V 2 + 2Ar 2 cos4°], +1 ]- 1 (20) 
is the 2D Greens function. The eigenstates of the latter operator within brackets are 



localized, with a localization length 27 1 Ro ~ A 2 //^. We therefore have G^(l,2) = 
(27r) _1 ln(_R / r i,2) at separations r 12 <C -Ro i n the weak-coupling limit, A ,A! — > oo (cf. 
ref. [7 1). A scale transformation f m — 1$ ■ x m of the 2n-dimensional integral above (fT9*j) 
yields the final result 

i^-KVAO^Ai/g 2 ]™ (21) 



for the asymptotic correlations of the superconducting order parameter across layers near 
zero temperature. The weakly coupled vortex-glass crosses over to a 3D vortex lattice 
threaded by lines of dislocations when the phase correlation length across layers, L^, exceeds 
the spacing between adjacent layers, d. This crossover occurs at a magnetic field 



B D (0) ~ (ya vx ) 2 (<£ /A 2 ; 



(22) 



near zero temperature, at which point the argument between brackets on the right-hand 
side of Eq. ()21|) is set to 1/e. The defective vortex lattice is decoupled across layers at 
perpendicular magnetic fields above B D (see Table EJ), where ^ < 

Consider again very weak Josephson coupling between adjacent layers, such that 1$ <ti 
Ai. Notice that this limit necessarily requires high perpendicular magnetic fields compared 
to the naive decoupling scale, ^o/Af, by the inequality a vx < 1$. Equation (}2*T]) then 
predicts short-range correlations of the superconducting order parameter across layers, with 
a correlation length that is less than the layer spacing d. Imagine next that the quenched 
disorder is reduced, such that 1$ 3> Ai. The argument in brackets on the right-hand side 
of Eq. (1211) then notably does not diverge towards positive infinity with the ratio Z^/Ai 
because of the logarithmic factor that originates from the 2D Greens function! Instead, it 
attains a maximum value of order unity at L ~ K\. Like in the cool-down from the vortex 
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liquid, Eq. these observations indicate that the correlation length for conventional 
superconducting order across layers does not diverge at perpendicular magnetic fields above 
the naive decoupling scale, B± > $ /Af. Unlike the case of vanishing thermal disorder 
(£0 — > 00) in Eq. (JSJ), however, the argument in brackets on the right-hand side of Eq. (J21j) 
diverges towards negative infinity with vanishing quenched disorder {1$ — > 00) because of 
the logarithmic factor. That divergence is spurious. The 2D Greens function (I20J) is given by 
G( 2 )(l,2) = (27r) _1 _K" ( r i,2/-Ro) m the limit — > 00, where cos^;°)/ +1 = 1. Here, K (x) is a 
modified Bessel function, and Rq = Ai/2 1 ^ 2 . Inspection of the original expression (|19j) for the 
autocorrelator across layers of the quenched superconducting order parameter then yields 
the asymptotic result lim^oo a„(i? /Ai) 2n = [(1 + [1 + (Ai/A ) 4 ] 1/2 )/4] n for that quantity 
as diverges. Notice that the latter argument raised to the power n instead saturates to a 
value that lies inside of the range [0.5,0.6], which is notably less than unity! No evidence 
for conventional superconducting order of the vortex lattice across a macroscopic number of 
layers therefore emerges from the above perturbative analysis to lowest non-trivial order in 
the Josephson coupling between layers, at B± > $ /A 2 . 



DISCUSSION AND CONCLUSIONS 



In conclusion, a duality analysis of the frustrated XY model for the mixed phase of layered 
superconductors with weak point defects finds that long-range vortex-glass order across 
layers emerges out of the vortex liquid at weak Josephson coupling between layers. This is 
consistent with recent Monte Carlo simulations of the same XY model that find evidence 
for a thermodynamic vortex glass phase jf| . It also potentially accounts for the recent 
observation of a thermodynamic vortex glass state in the mixed phase of high-temperature 
superconductors that show extreme layer anisotropy 9]. The analysis also indicates that 
the naive decoupling scale [25], $ /A 2 , serves as an upper bound for the stability of the 
Bragg glass phase as a function of perpendicular magnetic field in the extreme type-II limit. 
Previous theoretical work on layered superconductors predicts that the Bragg glass is stable 
to weak point pinning in general at the extreme type-II limit [11]. The discrepancy with the 
present work is likely due to the use there of a criterion for the destruction of the Bragg 
glass phase that is too stringent. In particular, the length along the field over which 
the vortex lattice tilts by a lattice constant is not divergent in ref. [11]. Also, the general 
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robustness of the Bragg glass predicted by ref. [ljj at weak pinning conflicts with the 
belief that the Bragg glass is generally unstable to invasion by dislocations in the limit of 
decoupled layers |13j|. Ai — > oo. A Bragg glass is also reported at fields beyond the naive 
decoupling scale in ref. [J, where the same XY model is studied numerically by Monte 
Carlo simulation. The discrepancy with the stability bound established here is likely due to 
a combination of finite-size effects and of intrinsic pinning by the grid in each 2D XY model 
(JIJ. The last effect has been neglected here throughout. Finally, Bragg peaks in neutron 
scattering that signal conventional vortex-lattice order at long range have been observed in 
the mixed phase of extremely layered high-temperature superconductors [12|, at fields below 
500 G. That threshold is consistent with the stability bound established here, 3>q/A^, if the 
Josephson penetration depth is bounded by A < 200 nm. Note that high layer anisotropy 
implies that the correction due to magnetic screening (A c ) suggested by ref. [25] can be 
ignored: Ai = A . 

The two theoretical results just reviewed depend critically on the existence of a vortex- 
glass state for isolated layers in the vicinity of zero temperature. Although recent experimen- 
tal determinations of the current-voltage characteristic in 2D arrays of Josephson junctions 
in weak magnetic field obtain evidence for melting of the 2D vortex lattice at transition 
temperatures T^ 2D ^ that are in fact much greater than the 2D melting temperature of the 



pristine vortex lattice |20|], theoretical arguments suggest that a perfectly conducting vortex 
glass can exist only at zero temperature in two dimensions [ll]]. Let us therefore consider 
the worst-case scenario, Tj 2D ) — > 0. The emergence of long-range vortex-glass order across 
layers from inside the weakly-coupled vortex liquid (fT2*|) survives this limit, since the 2D 
phase correlation length £ 2 r> remains divergent. Secondly, it is important to notice that the 
field equation (JT7j) used to obtain conventional phase correlations across layers (J2"T|) inside 
of the vortex glass is independent of the superfluid density pf ^ . This indicates that the 
stability bound in perpendicular magnetic field for the conventional vortex lattice, 3>o/Af, 
survives the limit Tf D ^ -> as well. 

The author thanks A. Koshelev for discussions. This work was performed under the 
auspices of the U.S. Department of Energy. 
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disorder index 


regime/phase 


(cos </>m+i) 


pi/Jz 


L^/d 




1 


Bragg Glass 


unity 


unity 


oo 


oo 


2 


Defective Vortex Lattice 


unity 


unity 


unity, or greater 


oo 


3 


Vortex Glass 


fraction 


fraction 


fraction 


oo 


4 


Critical Vortex Liquid 


fraction 





fraction 


unity, or greater 


5 


Decoupled Vortex Liquid 


fraction 





fraction 


fraction 



TABLE I: Listed are the conventional phase correlation length (L^) and the vortex-glass phase 
correlation length across equally spaced (d) layers, as well as the corresponding "cosine" and 
phase rigidity (see ref. for the various regimes found inside the mixed phase of an extremely 
type-II superconductor at weak Josephson coupling between layers, with weak point pinning. A 
horizontal line marks a true phase transition. 
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